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ABSTRACT 

We present Smoothed Particle Hydrodynamics (SPH) simulations of the response of 
gas discs to a spiral potential. These simulations show that the commonly observed 
spurs and feathering in spiral galaxies can be understood as being due to structures 
present in the spiral arms that are sheared by the divergent orbits in a spiral potential. 
Thus, dense molecular cloud-like structures generate the perpendicular spurs as they 
leave the spiral arms. Subsequent feathering occurs as spurs are further sheared into 
weaker parallel structures as they approach the next spiral passage. Self-gravity of the 
gas is not included in these simulations, stressing that these features are purely due to 
the hydrodynamics in spiral shocks. Instead, a necessary condition for this mechanism 
to work is that the gas need be relatively cold (1000 K or less) in order that the 
shock is sufficient to generate structure in the spiral arms, and such structure is not 
subsequently smoothed by the gas pressure. 

Key words: galaxies: spiral - hydrodynamics - galaxies: structure - galaxies: ISM 



1 INTRODUCTION 

Understanding the gas dynamics in spiral galaxies is an im- 
portant preconditi on if we are to compreh e nd the physics 
of star formation (lElmegreenI (I2002L [199^; iBonnell et all 
(|2006f) l and galaxy evolution. The origin of spurs and feath- 
ering is one phenomenon that has proved difficult to explain. 
These features are commonly observed in many spiral gjalax- 
ies, in cluding grand design galaxies such as M51 (jAalto et al.l 
(|l999h l. Spurs (particularly in M51) often show a correla- 
tion with HII regions, suggesting a link between these struc- 
tures and star forming giant molecular clouds. In fact the 
local star-forming region of Orion has been referred to as the 
Orion spur (Bokl959). Therefore, understanding the origin 
of spurs may offer insights into the star formation process. 

Spurs are seen as narrow dark lanes that extend per- 
pendicular to the spiral arms (|Lvndslll97Ci : IWeaveil 11970 '). 
They are particular ly striking in the HST images of M51 
(jScoville et aLlEooH) . especially towards the centre of the 
galaxy. The spurs extend over half the distance between 
the arms, and are often dotted with HII regions. Sur- 
veys indicate that spur s have pitch angles of 30 — 50"^ 
dRussell Roberts! Il992^ and widths comparable to those 
of spiral arms (lElmegreen|ll980h . Further ex amples of spurs 
are also shown in iBvrdI (jigS^ j (M31) and iKaufman et al.l 
(Il989) (M81). 

We denote feathering (iBalbusI Il988l: iKim OstrikeJ 
I2QQ2I ) by parallel filamentary features shifted from the spi- 
ral arms. These are often observed on the outskirts of spiral 
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galaxies (e.g. M74). Branches are more significant longer 
spurs often associated with resonances. Several simulations 
already have shown the bifurcation of spiral arms prim arily 
at th e 4:1 resonance of their potential (Patsis et al. 19941. 
ll997l : IChakrabarti et alll2003^ . 

Spurs have been suggested to form mainly through 
gravitational or magneto- Jeans instabilities. iBalbuJ (|l988f) 
demonstrated that transient gravitational stabilities develop 
preferentially nearly parallel and nearly perpendicular to the 
spiral arm. This introduces the appealing scenario where 
spurs and density perturbations along the arms are gener- 
ated by the same mechanism. Generally, dense regions along 
the course of the spiral arm will produce projections roughly 
perpendicu lar to the arm as a resul t of shear expansion. Sim- 
ulations bv lKim Ostrike3 (120021) show the formation and 
fragmentation of spurs in 2D magneto-hydrodynamical cal- 
culations of a shearing box. Dense clumps are produced in 
the arms principally through the magneto- Jeans instability, 
leading to the formation of spurs. 

Previous numerical analysis by iDwarkadas BalbuJ 

fl996) indicated that galactic flows through spiral arms are 
stable against purely hydrodynamic instabilities. However 
IWada Kodal ((20041) have looked again at the non- self- 
gravitating, non-magnetic case, performing global simula- 
tions of a spiral galaxy. They observe a rippling distortion of 
the shock front and the formation of spurs at larger pitch an- 
gles. These spurs are interpreted as a consequence of Kelvin- 
Helmholtz instabilities. 

We present hydrodynamic simulations which show the 
development of feathering and spurs in a spiral galaxy. 
These features are much more evident in our results than 
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previous purely hydrodynamical simulations. W e note that 
most previous stu dies, including IWada Kodal (|2004l) and 
towarkadas & Balbus (1996), assume a high sound speed. 
Our analysis reveals that the temperature of the disk has 
a crucial effect on the disk structure and the formation of 
spurs. We show that these features are due to the shearing of 
structure in the spiral arms. Furthermore, inhomogeneities 
transmitted between the arms lead to the formation of new 
structure in the next arm. 



2 CALCULATIONS 

We use the 3D smoothed partic le hydrody namics (SPH) 
code based on the version by Benz (iBenJ 199(1) . The smooth- 
ing length is allowed to vary with space and time, with the 
constraint that the typical number of neighbours for each 
particle is kept near N neigh ~ 50. Artificial viscosity is in- 
cluded with the standard parameters — 1 and (5 — 2 
dMonaghan Lattanzioll 19851 : lMonaghanlfl992r) . 

2.1 Flow through galactic potential 

The galactic potent ial includes a 4 armed spiral component, 
from lCox Gome d (^002). The symme tric components con- 
sists o f a logarithmic potential (e.g. iBinnev Tremaind 
(|l987f) ') that provides a flat rotation curve of vq — 
220 km s~^, and a p otential for the outer halo 
JCaldwell Ostrikeil[l98ir) . Parameters for the spiral part 
include the amplitude, 1 atom cm~^, pattern speed, 2 x 
10~^ rad yr~^ and pitch angle a = 15° which are compa- 
rable with the Milky Way. The pattern speed leads to a 
co-rotation radius of llkpc. 

Overall the disk is in equilibrium, as the rotational ve- 
locities of disk particles balance the centrifugal force from 
the potential. This paper only considers how hydrodynamic 
forces and galactic potential influence the flow. In particu- 
lar, self-gravity magnetic fields, heating, cooling or feedback 
from star formation are not included. 

2.2 Initial conditions 

Gas particles are initially chosen to occupy a region of ra- 
dius 5 kpc< r <10 kpc. The disk also has a scale height 
z ^ 100 pc. Particles are allocated positions and velocities 
determined from a 2D test particle run. The test particle run 
consists of 1 million particles initially distributed uniformly 
with circular velocities. They evolve for a couple of orbits 
subject to the galactic potential, to give a spiral density 
pattern with particles settled into their perturbed orbits. 
For the highest resolution run (4 million particles), 4.5 x 10^ 
test particles were split into 9 particles each. 

In the SPH initial conditions, the particles are given 
velocities in the z direction from a random Gaussian distri- 
bution of 2.5% of the orbital speed. The same magnitude ve- 
locity dispersion is also added to the x and y velocities. The 
gas is distributed uniformly on large scales with an average 
surface density of S 2 pc~^, although it is somewhat 
clumpy on smaller scales. The densities in the interarm re- 
gions range from 10~^^ g cm""^ to 10~^^ g cm""^. The total 
mass of the disk is 5 x 10^ M©. 

The number of particles is either 10^ or 4 x 10^. The 



higher resolution work was carried out on UKAFF (UK As- 
trophysical Fluids Facility). Simulations were run for be- 
tween 200 and 300 Myr, where the orbital period at 5 kpc 
is 150 Myr. No boundary conditions are applied to our cal- 
culations. 

All calculations are isothermal with temperatures of 50, 
10^, 10^, or 10^ K. Most previous simulati o ns and analysi s 
of spiral galaxies (e.g. Wada & Koda ( 2004*) ; "Gittins ("200j; 
iKim OstrikeJ ^2002); Dwarkadas & Balbus (1996)) as- 
sume a sound speed of approximately 10 km s~^ correspond- 
ing to the hotter (10^ K) component of the ISM. However, 
for GMCs to fo rm, gas ente ring the spiral arms must be cold 
atomic clouds rElmegreenl Eooal or pre-existing molecular 
gas (P ringle et al. 2001) . Since our initial aim of this work 
was to investigate molecular cloud formation (described in 
Dobbs et. al. 2006), lower ISM temperatures were adopted 
for these calculations. 



3 RESULTS 

The highest resolution run used 4 x 10^ particles, giving a 
resolution of 125 M© per particle. The temperature in this 
simulation was 50 K. Density plots at 4 different times are 
shown in Fig. 1. These and subsequent figures are shown in 
the rotating frame of the spiral potential. 

Initially there is a smooth spiral perturbation to an ap- 
proximately uniform disk. For the 50 K gas, the spiral arms 
are initially very thin. The density of the spiral arms in- 
creases with time and non-uniform structure develops along 
each arm. As this arm substructure grows, dense clumps 
shearing away from the spiral arms lead to interarm fea- 
tures. These features occur on the leading side of the spiral 
arm, leaving the trailing side smooth. Structure on the out- 
side of spiral arms is noticeable first in the inner parts of the 
disk (Fig. lb), and then spreads to larger radii. This is not 
surprising since the orbital period is less for smaller radii, so 
the comparative evolution of that part of the disk is faster. 
Both the arm and interarm regions become more disordered 
with time, and consequently the spiral arms become wider. 

Further runs have been performed at 100, 10"^ and 10^ 
K with 1 million particles and are shown in Fig. 2. As the 
temperature increases, the Mach number of the shock de- 
creases and much weaker spiral arms are produced. The gas 
is largely smooth at 10^ K due to the higher pressure, and 
the spiral arms exhibit no obvious substructure. 

The detailed structure of the spurs and feathering in 
the 50 K simulation is shown in Fig. 3. We see that the 
interarm structure increases with time from smaller scale 
structures to more pronounced larger scale spurs and feath- 
ering. Over time, spurs become more distinct with a regular 
spacing (?^ 0.7 kpc after 160 Myr in Fig. 3c). Furthermore, 
we see that there is a direct correlation between the spurs 
that leave perpendicular to the spiral arm and the feath- 
ering that approaches the next spiral arm. The feathering 
here has the appearance of closely spaced parallel lanes of 
gas with a much smaller pitch angle than apparent for spurs. 
The interarm structure in the 1000 K run is also shown in 
Fig. 2b. The spurs, and to a lesser extent the feathering, are 
still evident but are less pronounced. As noted above, the 
10^ K simulation does not display such interarm structure, 
but instead retains a smooth gas distribution (Fig. 2c). 
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Figure 2. Column density plots (g cm~^) when the temperature is a) 100, b) 10^ and c) 10^ K. The time corresponding for each figure 
is 220 Myr. Each plot shows a 5 kpc by 5 kpc section of the global simulation. The overall number of particles in these simulations is 1 
million. 



4 C. L. Dobbs and Ian A. Bonnell 



3.1 Formation of spurs 

In order to ascertain the origin of the interarm structure, we 
have traced the gas backwards in time from the spur that 
leaves the spiral arm to when the gas was located in the 
spiral arm itself. Figs. 4 and 5 show that there is a direct 
correlation between clumps in the spiral arms and the spurs 
that leave perpendicular to the arms. Thus we can conclude 
that the spurs are simply the shearing of clumps in the spi- 
ral arms, due to the divergent orbits of material that leave 
the arm. As structure (GMCs) is observed to exist in spiral 
arms, this provides a simple explanation for the presence of 
spurs. Further shearing then produces the feathering from 
the spurs. A similar scenario is discussed in Gittins (2004), 
who also describes the formation of spurs in SPH calcula- 
tions. 

The clumpy structure of the spiral arms is due to 
the dynamics as t he ISM passes through a spiral shock 
(iDobbs et al.ll2n0rf) . We interpret the growth of clumpy 
structure in the spiral arms in terms of the particles' change 
in angular momenta in the shock, which modifies the velocity 
phase space distribution of particles in the disk. The parti- 
cles enter the shock at non-uniform intervals of space or time 
(since the interarm regions are not homogeneous). As they 
subsequently gain and/or lose angular momentum in the 
spiral arms, the particles tend to group together in velocity 
space (Dobbs et al. 2006). Hence the inhomogeneities of the 
initial gas distribution are amplified in the spiral arms and 
clumps form. We estimate the spacing of these clumps from 
the time particles spend in the spiral arm (tarm ~ 6 x 10^ 
yea rs) and the typical velocity spread {v_\_ ^ 10 km s~^) 
dDobbs et al. 2006) to be L t ^ 600 pc. 

We plot the paths of five particles and their correspond- 
ing angular momenta in Fig. 6. Since the spiral arms are 
regions of significantly enhanced density, gas particles from 
a range of radii must interact in the shock. Particles enter- 
ing the shock on epicyclic paths are typically at the furthest 
extent of their orbit. They therefore encounter higher angu- 
lar momentum material already in the shock (e.g. red line. 
Fig. 6). This generally leads to a jump in the angular mo- 
mentum of the particles as they enter the shock (e.g. yellow 
line. Fig. 6). As particles travel along the spiral arm, they 
tend to lose angular momentum, since they are constrained 
by the spiral potential. It appears in Fig. 6 that particles 
entering the shock later (e.g. green line) show less change in 
angular momentum, possibly because they interact with gas 
in the shock which has similar angular momentum. Hence 
the path represented by the green line is less affected by the 
shock, and this particle spends less time in phase with the 
potential compared to the others. Eventually the gas leaves 
the spiral arm when it's angular momentum is too high to 
travel further radially inwards, and forms the spur shown in 
Figure 4.. 

This model also explains why less structure is apparent 
in the higher temperature simulations. At higher ISM tem- 
peratures, the shocks are weaker and they have less effect 
on the particles' angular momentum. The particles spend a 
shorter time in phase with the potential and their orbits are 
less perturbed. Overall the distribution of the angular mo- 
ment a of particles remains much more uniform (jPobbs et al.l 
I2OO6I I. Furthermore, the gas pressure can smooth out much 
of the structure during the crossing of the spiral arm. If the 



structure in the spiral arms is removed, then there is noth- 
ing to shear into the subsequent spurs and feathering in the 
interarm region. For example, the 10^ K run has a sound 
speed of ~ 10 km s~^ which is comparable to the velocity 
of the gas perpendicular to the spiral arm (in the rotating 
spiral potential). Thus, structures of size scales ^ the width 
of the arm can be smoothed by the internal pressure before 
leaving the spiral arm. 

An alternative explanation is that the spiral arm 
structure arises fr om Kelvin Helmholtz instabilities 
(IWada Kodal l2004ri . We find that the Reynolds numbers 
for our simulations are relatively low {Re ~ 500 for 100 K 
gas due to the numerical resolution) implying that the res- 
olution is insufficie nt for Kelvin H e lmhol tz instabilities to 
occur. Interestingly, IWada Kod3 (l2004l) find spurs form- 
ing from hot (10^) gas, although their analysis likewise shows 
that increasing the Mach number of the shock (correspond- 
ing to a lower temperature) increases the susceptibility of 
gas disks to hydrodynamic instabilities. 



4 CONCLUSION 

We have performed numerical simulations of disks in spi- 
ral galaxies which demonstrate the widespread formation of 
spurs and feathering. The overall appearance of the galactic 
disk, enhanced by these features, is remarkably similar to 
observations of galaxies such as M51 and M81. In addition, 
the properties of individual spurs (pitch angle, dimensions) 
are comparable with observations. 

The formation of spurs and feathering in the interarm 
regions of spiral galaxies is shown to be linked to the pres- 
ence of structure in the spiral arm s. These findings concur 
with similar work by Gittins ("2004*) . Upon leaving the spiral 
arms, dense concentrations of gas are sheared by diverging 
orbits. Spurs form perpendicular to the spiral arms, which 
later become feathering that impacts into the next spiral 
arm. The formation of spurs and GMCs are thus directly 
related, since these denser structures in the spiral arms are 
the sites of molecular clouds. 

The necessary conditions for interarm features to form 
are density inhomogeneities and relatively cold gas (T < 
1000 K). A clumpy ISM allows structure growth in cold spi- 
ral shocks. Such structure formation in the spiral arms, and 
therefore spur formation, does not occur in our simulations 
when the gas is hot. 
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Figure 5. Formation of a spur (right hand plot, after 220 Myr) from a dense clump of gas in the spiral arm (left hand plot, after 180 
Myr). The temperature is 1000 K and both plots are 2 kpc by 2 kpc. 
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Figure 6. The orbits of 5 selected particles are displayed pre- 
and post-shock (top) where the xy-coordinates correspond to a 
fixed Cartesian grid centred at the midpoint of the disk. The 
dotted line indicates the position of the shock front, and the orbits 
are taken from the rest frame of the potential. Boxes indicate 
the location of the five particles at three corresponding times. 
The final position of the particles coincide with the spur shown 
in Figure 4 (right). The corresponding angular momenta of the 
particles are also plotted (bottom). Particles are travelling from 
right to left in both plots. 
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